Realization of Bose-Einstein condensation of dilute magnons in TlCuCla 
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The recent observation [Oosawa et al. J. Phys. : Condens. 
Matter 11, 265 (1999)] of the field- induced Neel ordering in 
a spin-gap magnetic compound TlCuCls is interpreted as a 
Bose-Einstein Condensation of magnons. A mean-field cal- 
culation based on this picture is shown to describe well the 
temperature dependence of the magnetization. The present 
system opens a new area for studying Bose-Einstein conden- 
sation of thermodynamically large number of particles in a 
grand-canonical ensemble. 

PACS number 75.10.Jm 



The Bose-Einstein Condensation (BEC) is one of the 
most exotic phenomena predicted by quantum mechan- 
ics Although the superfluid transition of Helium 4 
may be regarded as a BEC, it is influenced strongly by 
the interaction and is much different from condensation 
of ideal or dilute Bose gas. Recently, there has been a 
renewed interest in BEC, because the realization of BEC 
by ultracooling of dilute atoms has become possible [Q. 
While the BEC of ultracooled atoms is of great interest, 
there are various experimental limitations. On the other 
hand, it has been known for a long time that a quantum 
spin system can be mapped to an interacting Bose gas, 
and that the off-diagonal long-range order which charac- 
terizes BEC corresponds to a long-range magnetic order 
in the spin system ||^. It is then possible to tune the 
density of bosons (magnons) by a magnetic field to ob- 
serve BEC of dilute bosons. However, such an attempt 
has been apparently lacking. 

In this letter, we argue that BEC of dilute bosons in 
a thermodynamic number ^ 10^° is realized in a recent 
high-field experiment on TlCuCla B, which is composed 
of chemical double chain of CU2CI6 [|^,^. The compound 
has an excitation gap A/fce ~ 7.5K above the siiiglet 
ground state, in the absence of the magnetic field |@,||,@|. 
The origin of the gap may be attributed to the antiferro- 
magnetic dimer coupling in the double chain. When the 
external field Hg = A/{gfj,B) to the gap is apphed, the 
gap collapses. At finite temperature, the "collapse" of the 
gap at Hg does not give a singularity because thermal ex- 
citations exist even if iJ < Hg. However, there seems to 
be a phase transition due to the interchain interactions at 
higher field H = He > Hg , which depends on the temper- 
ature. In Ref. [H the phase transition was identified as a 
long-range magnetic ordering, and was compared with a 
mean-field theory (MET) based on a dimer model. 
While the dimer MET does predict the field-induced or- 



dering, the experimental features were not well repro- 
duced. In particular, it predicts almost flat dependence 
of the critical temperature Tc on the magnetic field, while 
in the experiment Tc depends on the magnetic field by 
a power law T;!' H - Hg (see Eig. 1). Moreover, it 
predicts almost constant magnetization for T < Tc and 
concave magnetization for T > Tc, as a function of tem- 
perature T . However, in the experiment, magnetization 
was found to increase as decreasing T below Tc and it is 
a convex function of T for T > Tc (see Fig. 2). 

We will show that the transition is rather well de- 
scribed as BEC of magnons. While the details of ex- 
change interactions in TlCuCla are not known yet, exci- 
tations above the singlet ground state generally can be 
treated as a collection of bosonic particles - magnons . 
If the exchange interaction is isotropic, which seems to 
be the case in TlCuCla , the number of magnons are con- 
served in a short timescale (but not conserved in a longer 
timescale.) We assume that magnons carry spin 1, as 
generally expected. 

Under a magnetic field H ~ Hg, the magnons with 
S'^ — 1 can be created by small energy. Thus, at low 
temperatures T <C A and H ^ Hg , we can consider only 
those magnons. The chemical potential of the magnons 
are given by /i = gfi^iH — Hg). The total number of 
magnons N is associated with the total magnetization 
M through M = gii-^N . If the magnons are free bosons, 
the number of magnons would be infinite for H > Hg. 
However, in the spin system, magnons cannot occupy 
the same sites and thus there is a hard-core type inter- 
action between them. The interaction keeps the number 
of magnons to be finite. 

The transverse components of the exchange interac- 
tions give rise to hopping of the magnons, while the longi- 
tudinal component give rise to the interaction. Although 
the exchange interaction and thus hopping might be com- 
plicated, generically the dispersion relation of a magnon 
is quadratic near the bottom. Thus the low-energy ef- 
fective Hamiltonian for the (S^ ~ 1) magnons are given 
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Here the momentum k is measured from the minimum of 
the magnon dispersion. Eor simplicity, we do not consider 
the case where the magnon dispersion has more than one 
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minimum . The effective masses nia is related to the 
curvature of the dispersion relation in the direction of a. 
By a rescaling of momentum, we may consider isotropic 
effective Hamiltonian instead. This is nothing but the 
non-relativistic bosons with a short-range interaction. 

Moreover, in the low-density and low-temperature 
limit, only the two-particle interaction is important and 
it can be replaced by delta-function interaction v{q) ^ vq. 
Thus the effective Hamiltonian is given by 
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This effective Hamiltonian can be derived from some spe- 
cific models However, we emphasize that it is 
universal at the low-temperature and low magnon den- 
sity limit, and does not depend on details of the exchange 
interaction. 

Since the number of magnon is actually not conserved 
due to the small effects neglected in the Hamiltonian, 
we have a grand canonical ensemble of the bosons. The 
"chemical potential" can be controlled precisely by tun- 
ing the magnetic field. When the chemical potential be- 
comes larger than a critical value, the system undergoes a 
BEC. Thus the spin-gap system in general would provide 
a great opportunity to study BEC in a grand canonical 
ensemble, with thermodynamically large number of par- 
ticles. 

The idea that BEC is induced by the magnetic field in 
a spin-gap system has appeared several times. There was 
a discussion of (quasi-) Bose condensation in a Haldane 
gap system under a magnetic field fl^ . However, there 
is no BEC at finite temperature in a one-dimensional 
system. On the other hand, the experiments on Hal- 
dane gap systems are often affected by the anisotropy 
and the staggered g-tensor, which wipe out the BEC. 
Giamarchi and Tsvelik have recently discussed the 
three-dimensional ordering in coupled ladders in connec- 
tion with BEC. However, as far as we know, there has 
been no experimental observation of the magnon BEC 
induced by an applied field. 

We first consider the normal (non-condensed) phase. 
Within the Hartree-Fock (HF) approximation, the mo- 
mentum distribution of the magnons is given by |15| 
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with Ek = h k /2m and fics = 1^ ^ 2von. The magnon 
density n = N/Nd {Nd is the total number of the dimer 
pairs) has to be determined self-consistently by 
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where z = e'^'''^" is the fugacity, A = {27^1%^ /mk-BTf''^ 
is the thermal de Broglie wavelength, and gn{z) = 
Sz^i 2'/^" is the Bose-Einstein function. BEC occurs 
when the effective chemical potential fj,cs vanishes so that 
fi = 2von. Setting z = 1 in (4) gives the temperature de- 
pendence of the critical value of the chemical potential 
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This implies that the temperature dependence of the crit- 
ical magnetic field at low temperatures is Hc{T) — Hg cx 
T^/^. This power-law dependence is independent of the 
interaction parameter uq. 

When /i > /ic, one has the macroscopic condensate or- 
der parameter (ao) = VlVce*^ ^ 0, where Nc is the total 
number of the condensate magnons. In terms of the orig- 
inal spin system, this means that there is a (staggered) 
transverse magnetization component m^ — 5/iB \/nc/2 
with Tic = Nc/Nd- Within the Hartree-Fock-Popov 
(HFP) approximation, the condensate density is deter- 
mined by |l6[ 



fJ. = vqUc + 2vQh, 
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where n = n ~ Uc is the density of the non-condensed 
magnons, which is given by 
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where we have used the HFP energy spectrum E^ — 
\Js\^ 2ekVonc and the Bose distribution fB{Ek) — 
l/{e^^'' — 1). The first term of (7) represents the 
depletion of the condensate due to interaction be- 
tween magnons, which reduces to the ground-state non- 
condensate density at T — » 0. The second term is 
the contribution from thermally excited non-condensate 
magnons, which vanishes at T ^ 0. Eq.(6) is to be solved 
self-consistently in conjunction with Eq.(7). Then the to- 
tal magnon density is given by 



n — Tic + n — n. 
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In particular, the magnon density at T — > is given by 
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If we ignore the deviation of z from 1 for T ^ Tc, we 
obtain a simple result (using vqUc = in (7) for T < 
Tc) pi: 
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Thus it predicts the cusp-Uke minimum of the magnon 
density (magnetization) at T — Tc. In contrast, the 
dimer MFT js) predicts a constant magnetization below 
Tc. 

Figure ^ shows the observed low-temperature magne- 
tizations of TlCuCla at various external fields for H \\ b. 
We can see the cusp-like anomaly at the transition tem- 
perature, as predicted by the present theory. The simi- 
lar temperature dependence of the magnetization can be 
observed for _L (1,0,2) Thus the main qualitative 
feature of the temperature dependence of the magnetiza- 
tion, which cannot be understood in the dimer MFT, is 
captured by the magnon BEC picture. The increase of n 
for decreasing T below Tc is due to condensation of the 
bosons; the cusp shape of the magnetization curve ob- 
served in the experiment can be regarded as an evidence 
of the magnon BEC. We note that, in the range of the 
experiment, the magnon density is of order of 10"'^ and 
is consistent with the assumption of diluteness. 

However, the approximation ( [l0| ) does not precisely re- 
produce the experimental result. In particular, it predicts 
independence of n on the applied field fi for T > Tc while 
the dependence was observed experimentally. Part of the 
discrepancy may be due to the approximation z = I. 
Actually, even in the HE framework, the approximation 
z = 1 cannot be justified. In Fig. 3, we plot the temper- 
ature dependence of the total density n above and below 
the transition temperature Tc obtained by solving the 
self-consistency equations (4) and (8) numerically. The 
interaction parameter vq and the effective mass m are 
estimated from the experimental data as wq/^b ~ 400K 
and mk^/Ti'^ w 0.025K^^. The self-consistent calcula- 
tion does predict the total density n dependent on the 
applied field for T > Tc, which is qualitatively consistent 
with the experiment. In Fig. |4| we also plot the temper- 
ature dependence of the staggered transverse magnetiza- 
tion component m_L. Direct measurements of m± using 
neutron diffraction are in fact intended. 

We see a discontinuity in magnon density (magnetiza- 
tion) at the transition point. This is because our HEP 
approximation is inappropriate in the critical region, and 
leads to an unphysical jump in the condensate density Hc 
(for detailed discussion, see jl^). In the vicinity of the 
critical point, the HEP approximation eventually breaks 
down; the critical behavior then belongs to the so-called 
3D XY universality class [^. On this ground, in the 
vicinity of Tc, the transverse magnetization m± is ex- 
pected as m± cx {Tc — T)^ , where (3 ~ 0.35. 

Figure ^ shows the experimentally determined mag- 
netic phase diagram of TlCuCls. We fit the phase bound- 
ary He as a function as a temperature Tc with the fol- 
lowing formula: 



The best fitting is obtained with {g/2)Hc{0) = 5.61T 
and = 2.2 [|l9|. The obtained exponent (j> = 2.2 dis- 
agrees somewhat with the HE approximation (|^) which 
gives (j) = 3/2. We note that, exactly z = 1 holds at 
the transition point, and thus ^ = 3/2 is a definite con- 
clusion within the HE framework. On the other hand, 
the dimer MFT predicts Hc{T) to be exponentially flat 
at low temperature Q . The observed power-law depen- 
dence is qualitatively consistent with the magnon BEC 
picture, compared to the dimer MFT. 

As discussed above, our mean-field analysis for a dilute 
Bose gas is not reliable in the critical region, and thus 
the discrepancies with the experiment may be attributed 
to the fluctuation effects. More precise description of 
the experiment near the critical point therefore requires 
the inclusion of the fluctuation effects. Furthermore, in 
the experiment there may be other effects that were ig- 
nored in the effective Hamiltonian such as impurities. 
These will be interesting problems to be studied in the 
future. 

To conclude, we believe that the essential feature of the 
experimental observation on TlCuCla, which cannot be 
understood in the traditional dimer MET, is captured by 
the magnon BEC picture. The present system provides 
the first clear experimental observation of field-induced 
magnon BEC, with thermodynamically large number of 
particles. It opens a new area of BEC research in grand 
canonical ensemble with an easily tunable chemical po- 
tential (magnetic field). Similar BEC of magnons would 
be observed in other magnetic materials in the vicinity 
of the gapped phase, which may be the singlet ground 
state due to large single-ion anisotropy [Q, the com- 
pletely polarized state (|jl^, or the "plateau" phase in 
the middle of the magnetization curve j2l] . An essential 
requirement for observing BEC is that the system has 
the rotationally invariance about the direction of the ap- 
plied magnetic field, so that the number of magnons is 
(approximately) conserved. 
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FIG. 2. The low-temperature magnetizations of TlCuCls 
measured at various external fields for 11 b. 
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FIG. 3. The temperature dependence of the magnetiza- 
tion. We have used vo/ks = 400K, mkn/Ti'^ = 0.025K"^ and 
g = 2.06. 
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FIG. 1. The phase diagram in TlCuCls. The soUd line 
denotes the fitting with Eq.(ll) using {g/2)Hc{0) = 5.61T 
and (b = 2.2. 



4 




5 



